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Abstract 

The extension of the minimal standard model by three right-handed sterile neutrinos with 
masses smaller than the electroweak scale (i^MSM) is discussed in a Qq flavor symmetry frame- 
work. The lightness of the keV sterile neutrino and the near mass degeneracy of two heavier 
sterile neutrinos are naturally explained by exploiting group properties of Qe- A normal hi- 
erarchical mass spectrum and an approximately fi-T symmetric mass matrix are predicted for 
three active neutrinos. Nonzero ^13 can be obtained together with a deviation of 823 from the 



maximality, where both mixing angles are consistent with the latest global data including T2K 
and MINOS results. Furthermore, the tiny active-sterile mixing is related to the mass ratio 
between the lightest active and lightest sterile neutrinos. 
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I. INTRODUCTION 



Both the estabhshment of neutrino oscillation phenomena and the evidence of nonlu- 
minous dark matter (DM) demand physics beyond the standard model (SM). On the one 
hand, compelling evidences from current solar, atmospheric, reactor and accelerator neu- 
trino experiments have told us that neutrinos are massive and lepton flavors are mixed 
On the other hand, various cosmological observations have revealed that DM is five times 
more abundant than normal matter and accounts for about one quarter of the Universe 



(i.e., n^M ^ 0.227 ±0.014 



2j). To explain both problems above, an extension of the SM 



by three right-handed sterile neutrinos with masses smaller than the electroweak scale 



(z/MSM) was first proposed by Asaka et al. [3j. The smallness of active neutrino masses 
is described by the canonical seesaw mechanism [4] and one light right-handed sterile neu- 
trino at the keV scale acts as a candidate of warm dark matter (WDM). Moreover^the 

tin 

model can explain the baryon asymmetry in the Universe through the oscillations [5|, 16| 
of two heavier right-handed sterile neutrinos with masses at the GeV scale. 

Despite the above phenomenological successes, one unsatisfactory point in the z/MSM is 
the lack of a natural explanation for large mass splitting between the keV sterile neutrino 
and the heavier ones. Moreover, the oscillation mechanism for the baryon asymmetry 
demands strong mass degeneracy between the heavier sterile neutrinos, but the i/MSM is 
impotent to be able to explain its origin. In order to resolve theses issues, some interesting 
ideas have been proposed in the context of a U(T) flavor symmetry the Froggatt-Nielsen 



mechanism 



□ 

the split seesaw mechanism [9|] and grand unified theories [10|]. In the 
present work, we introduce a non-Abelian discrete flavor symmetry and try to understand 
the mass splitting and degeneracy due to group properties of the flavor symmetry. The 
mass degeneracy of two heaver sterile neutrinos could be interpreted as a sign that they 
constitute a doublet representation of the flavor symmetry, while it may be natural to 
assign a singlet representation to the keV sterile neutrino. Furthermore, if the singlet 
representation is a complex one, we can prohibit a bare mass term of the keV sterile 
neutrino because of the Majorana nature and may be able to generate a suppressed mass 
term from higher-dimensional operators. Inspired by these clues, we employ the Qq group 
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TABLE I: Particle content and charge assignments. 

as our flavor symmetry^ since it is the smallest finite group^ which contains both a complex 
singlet and a real doublet representations. In addition to Qq, we introduce two auxiliary 
Zjv symmetries in order to handle the order of magnitude of some small parameters as well 
as to forbid unwanted terms. We find that our model predicts a normal hierarchical mass 
spectrum and an approximately fi-r symmetric mass matrix for three active neutrinos. 

The remaining part of this paper is organized as follows. In section II we present the 
framework of the i^MSM with a. Qq x Z3 x Z2 fiavor symmetry. Section III is devoted to 
the realization of the seesaw mechanism and section IV to the diagonalization and the 
resulting neutrino masses and mixing matrix. A numerical analysis with the focus on 
nonzero ^13 is also illustrated in section IV. Finally, we give a summary in section V. 



II. THE QqX Z3X Z2 MODEL 

We introduce three right-handed sterile neutrinos and gauge-singlet fiavon fields 

Sx, Sy, Sz and D with a. Qq x Z3 x Z2 fiavor symmetry. We assign a Qe-singlet (- 
doublet) for the first (second and third) generation of fermions, and the SM Higgs is 
assumed to be invariant under all the fiavor symmetries. The particle content and charge 
assignments for each symmetry are summarized in Table I, and the basic group theory of 
Qq is reviewed in the appendix. Because of the symmetries, there are no renormalizable 
Yukawa interactions in the charged lepton sector, and the charged lepton masses follow 
from the higher dimensional operators: 

Ce = ^(L^HE^\,„S, + ^iL^HE^)^„Sy + ^ L,HE,Sl + h.c. , (1) 



1 



Previous ati 



dies about Qq can be found in Ref. [UJ 



^ See Ref. 12 1 for the classification of discrete groups upto g = 31, and Ref. for the recent 



reviews 



of non-Abchan discrete flavor symmetries. 



where we have specified the Lagrangian upto the next-to-leading order level, and the 
charged lepton mass matrix is diagonal upto this order. The subscripts beside (■ ■ ■) in- 
dicate the required representations to make the terms invariant under Qq. In the 
neutrino sector, Ni is assumed to be a complex representation of Qe? while the other 
leptons are real representations. Consequently, only Sx and S'y, which are complex repre- 
sentations, can reproduce a Dirac and Majorana mass terms for A'^i, and those interactions 
are suppressed due to the Z^r symmetries: 



+ — L,HN,Sl + ^ (L^HN^)^,S.,S;S, + h.c. 

TTlh TTl 



h.c. . 



(2) 
(3) 



Note that in the above expressions we have omitted some terms whose contributions can 
be embedded into others and implicitly assumed a mechanism which supplies ma^b,c = 
(9(1) GeV, e.g., the spontaneous breaking of a lepton number symmetry at a GeV 
scale. 

We define the vacuum expectation values(VEVs) of neutral scalars as 



i/O) = ^; = 174 GeV, {S,) = s,, (Sy) = Sy, (5,) = s„ (D) = {d,,d2), 



(4) 



leading to the charged lepton masses 

2 1 



i 1 

^ j YeV, ^ti = O^xSx + YySy) V, = O^xSx " YySy) V 



the right-handed Majorana neutrino mass matrix 
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and the Dirac mass matrix, 

^ ad2 adi ^ 



1^ nibdl J 



/3di 'ys, 
\ (3d2 J 



I Ssl 



A3 








aSxSyS^ 



^SxSySz 



V. 



/ 



(5) 



(6) 



(7) 
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Notice that, in Eq. we have the same fine-tuning problem as that in for obtaining 
rrin ^ mr. Nevertheless, we will not tackle this problem and only focus on the neutrino 
sector in what follows. We also note that although the second terms of Mji and are 
strongly suppressed by l/A^ compared with the first terms, we keep mc, rrib and 5 in our 
discussions because they will be important when we discuss the sterile neutrino masses 
and active-sterile mixing. In contrast, e in M£, only contributes to the masses and mixing 
of three active neutrinos, and its effects are negligibly small in comparison with the first 
term. Thus, we shall ignore e. 



III. SEESAW MECHANISM 



Let us move on to the diagonal basis of M^. The diagonalization can approximately 
be done by the 45° rotation in the 2-3 plane, and three sterile neutrino masses are found 
to be 
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Ml is suppressed with l/A^, and thus it is the candidate of WDM, while M2 and M3 
are nearly degenerate. Interestingly, the order of Mi and that of the mass difference 
between M2 and M3 are the same, which may turn out to be a key ingredient when 
one considers the baryon asymmetry in the Universe 6|]. Nevertheless, we shall naively 
assume M2 = M3 = and do not consider the baryon asymmetry in what follows since 
detailed studies of the baryon asymmetry go beyond the scope of this paper. Because of 
Ml <^ M2,3, the masses of the light neutrinos are obtained by integrating out only M2 
and M3, yielding the following 4x4 effective mass matrix: 

/ 



m; 



4x4 




(9) 
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(10) 



a'{2(3' - 7') 2(3'i 
y -a'{2/3' - 7') -2/3'y 
will later end up the mass matrix of three active neutrinos, and A = {6'v, 0, 0)"^ stands 
for the mixing effect of active-sterile neutrinos. In the above expressions, we have defined 
di + d2 = 2d and di — d2 = 2edd and embedded the suppression factor 1/A into the 
couplings such that 



6' = 6 



d 



a 



a- 



/3' = /3t, 7^ 
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a; ' A ' " "A ' ' 'A 

As one can see, Mj^^^ is fi-r symmetric if = 0{di = ^2) holds. 

Here, let us roughly estimate the magnitude of model parameters. From the charged 
lepton sector, we obtain s^^y/A ^ (2 • ■ - 5) x 10"'^ for Ye^x,y = Suppose s^/A = 10"'^ 

and Sx — Sy c:^ Sz — di c:^ d2 for simplicity, then = 0{1) GeV results in Mi = 0{1) keV 
for the lightest sterile neutrino. Since M2,3 ~ rria and = 0{1) GeV, {a', /?', 7') need 
to be (9(10"^'^) in order to reproduce realistic active neutrino masses = (9(10~^) eV, 
and they correspond to {a, /3, 7) = (9(10~^'^). If we assume the same value for 6 as well, 
then we gain 6' = O{10-^^-^). 

In the approximation of 6'v <^ Mi, the mass matrix M^^^ in Eq. 
diagonalized by a 4 x 4 neutrino mixing matrix parametrized as 

I \ I ± — 111 I 11 1 ll^AU 

ViVo 



can further be 





(12) 



where i? is a 3 x 1 mixing matrix which induces the active-sterile mixing angles 

R ^ A/M^ = 0, 0)^ , (13) 

and Va is a unitary 3x3 mixing matrix diagonalizing the mass matrix of three active 
neutrinos. We can use Vi to eliminate the A and A"^ terms of M^^'^, 

Mf.'^^ + RMiR^ 
Ml 

and Va can be determined by the relation Va(M^^^ + RM^R^)Vl = diag {A^, Ag, A3}. 
Because the order of RM^R^ [i.e., ^'^f^/M^ = O{10^^) eV] is much smaller than that of 
M^^^, we can safely neglect its effects in the active neutrino part. The only exception is 
the contribution to A^ because it is vanishing when M^^^ is taken into account alone. 



(14) 
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IV. NEUTRINO MASSES AND MIXING 



ACTIVE NEUTRINO MASSES AND MIXING 



As we mentioned, if = holds, M^^^ is /i-r symmetric and it results in the vanishing 

act, these predictions are roughly compatible 
-|l7l|. which indicate a small 6*^3 and a nearly 



and maximal ^23 for active neutrinos. 



n 



15 



with the recent neutrino oscillation data 
maximal 6'23. Thus, Ed may be expected to be not so large, so that we here treat 
as a small parameter and employ the perturbation calculations. The leading (i.e., /i-r 
symmetric) term can be diagonalized by 
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where cq = cos 6 and sq = sinO, and the mixing parameters are given by 
p = arga', a = arg(/3' + 7'), sin6' v^lo' I 



y^|^' + 7'|2 + 2|a'| 



(16) 



together with the eigenvalues 



A? = 0, A 



|/3' + 7f + 2|a^ 



J\2 



(/3'-7T . 



(17) 



Notice that p and a are unphysical phases and do not affect any observables, though they 
appear in the following expressions. After including the correction term and doing per- 
turbative calculations, we get the neutrino mixing angles in the standard parametrization 

as 



sin 6^12 - sin6' + 0{el) 
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and the Jarlskog invariant parameter J as 

y2|a'(2/3'-7')l 



(19) 



where (f) = arg(2a' + /3' + 27'). From Eq. fll8l) . one can check the recovery of the ^-t 
symmetry (^^3 = 0° and ^23 = 45°) in the hmit of Ed = 0. The corrections to the three 
eigenvalues are vanishing in the order of 0{e^), so the active neutrino masses are obtained 
by taking absolute values for Eq. (1171) . and this model predicts a normal mass hierarchy 
for three active neutrinos. Note that the vanishing of Ai should be kept in all orders of 
perturbations and is a generic property of the minimal seesaw models 18|]. A nonvanishing 



mass of Ui can be generated from the lightest sterile neutrino contribution: 

Ai ^ 6'\yM,, (20) 

but the corresponding effects are negligibly small and can be safely ignored for all the 
other mass and mixing parameters. 

B. NUMERICAL ANALYSIS 

Instead of perturbative calculations, we here numerically diagonalize Eq. (fTOl) and 



compute 6'i3, 623 and J. From the recent global analysis 17|] of the neutrino oscillation 



data, in our calculations, we refer to the following best-fit values and 3a error bounds: 

Ami, = (7.59^°:^°) X 10"^ eV^ Ami = (2.50l°:i) x 10"^ eV^ 

sin^ = 0.312l°:°^«, sin^ = 0.52lO:J2, sin^ = 0M3t',f,l (21) 

for the normal neutrino mass hierarchy. In Figure [H we plot sin^ ^^3 as functions of 
sin^ 6'23 (left panel) and J (right panel) with the 3a constraints of Am"^,, Am"^, and 
sin^ 9,2- Besides, the 3a bound of sin^ 6*23 is also imposed in the sin^ 9,^ — J plane (right 
panel). As one can see from the figure, the predicted regions can be within the 3a ranges. 



and sin 6*^3 can deviate from 0, which is favored by the recent T2K [l9|] and MINOS 20 1 



results. However, a large sin^ 9,^ is always accompanied with a large deviation of sin^ ^23 
from 0.5. For instance, the best-fit value of sin^ ^^^3 can be accounted for at around 
I sin^ 6*23 — 0.5| ~ 0.11, but it is almost the edge of the 3a bound^. Notice that we have 
checked that the h-t symmetry breaking parameter e^, which is defined below Eq. (fTOl) . 



^ We refer to Ref. 
maximality. 



2ll | for the realizations of a large 9,^ together with a small deviation of from the 



8 



0.025 




0.025 



0.65 




0.03 



FIG. 1: sin ^2.3 ^23 panel) and sin 9^^ vs tJ (right panel), where the horizontal dotted line 

displays the best-fit value of sin^ while the vertical dashed lines express the 3a upper and lower bounds 
of sin2 6'23 from Eq. (^1]) . 



is at most 0.4 and that the analytical expressions in Eqs. (H 
agree with the numerical results in Figure [H 



and (flSll approximately 



C. ACTIVE-STERILE MIXING 



From Eq. (fT3|) and the parameter estimates in section III, the active-sterile mixing 
angles Gj defined by the elements of the R matrix (i.e., 0j ^ l-Rjil) '^^^ be derived as 

Ai 



Ml 



10 



-11 



(22) 



and 62 ^ G| ^ 0, which are well below the upper bounds from astrophysical and cosmo- 



22|. Furthermore, it is also consistent with the requirement of correct 



logical observations 

DM abundance for the mechanism of resonant active-sterile oscillations with nonzero lep- 



ton asymmetries 
anism 



231]. To achieve the right DM abundance with the nonresonant mech- 



24i\ . we need 6 to be one order of magnitude larger than (a, /3, 7) in order to 



achieve 0^ ^ 10 ^ [223. These tiny active-sterile mixing angles make the detection of 



the WDM particle rather dim and remote with both the X-ray observations 



and the 



captures on beta-decaying or electron-capture-decaying nuclei j26l. |27|. 
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V. CONCLUSIONS 



In this work, we have proposed a flavor symmetry reahzation for the z/MSM in 
the presence of two auxihary symmetries and succeeded in naturally explaining the 
lightness of the keV sterile neutrino and the mass degeneracy of the two heavier sterile 
neutrinos. A normal hierarchical mass spectrum and an approximately /i-r symmetric 
mass matrix are predicted for three active neutrinos. Nonzero ^13 can be obtained together 
with a deviation of 6*23 from the maximality, where both mixing angles are consistent with 
the latest global data including T2K and MINOS results. Finally, we have derived a tiny 
active-sterile mixing related to the mass ratio between the lightest active and lightest 
sterile neutrinos. 

The i/MSM can explain the active neutrino masses, the candidate of dark matter and 
the baryon asymmetry in the Universe in a unifled and elegant way. Although there are 
already some models in which the z/MSM is extended by flavor symmetries, our model has 
more direct connections with the masses and mixing patterns of three active neutrinos. 



Our realization can also be modifled to accommodate the eV scale sterile neutrinos 



28 



which are more or less hinted at by current experimental 29|] and cosmological 30| data. 



We shall examine this case with a speciflc flavor model elsewhere 



3l|. 
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Appendix A: Basics of 

Qq consists of four singlet and two doublet irreducible representations, 

I n III I , . \ 

1, 1 , 1 , 1 , 2, 2 , (Al) 
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TABLE II: A table of the tensor products of Qq. 



and twelve elements, 



Rq, i?g, ■ ■ ■ i?g, Pq, RqPq, RIPq, ■ ■ ■ RIPq, 



(A2) 



where E stands for the unit matrix. The representation matrices of Rq and Pq for each 
representation are give by 
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with ojg = exp {^^}- Note that 1, 1 and 2' are real representations, 2 is a pseudoreal 
representation and 1 " = (l")* are complex representations. The tensor products among 
the irreducible representations are summarized in Table II. Especially, the products of 



11 



two doublets are defined as follows. 
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